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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 

Two  mathematical  descriptions  of  the  strain-rate  dependence  of  materials  response 
to  suddenly  applied  one-dimensional  strain  loads  were  incorporated  into  the  PUFF 
wave  propagation  computer  codes.  The  models  are  quite  different  in  that  one 
uses  a  macroscopic  relaxation  time  to  govern  the  stress-relaxation  rates,  while 
the  other  uses  the  microscopic  theories  of  dislocation  dynamics.  The  results  of 
several  sample  calculations  are  discussed,  and  the  changes  required  in  the 
programming  sequence  are  given.  Some  comments  are  provided  concerning  the 
validity  and  usefulness  of  these  models,  and  a  brief  analysis  of  the  effects  of 
parametric  variations  is  presented. 
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SECTION  I 
INTRODUCTION 

The  PUFF  66  and  P  PUFF  66  computer  codes  were  developed  to  describe  the 
shocks  resulting  from  thermally  and  mechanically  Induced  loads  on  semi-infinite 
slabs  of  material  (Reference  1).  These  codes,  which  compute  one-dimensional 
motion,  are  based  upon  the  principles  of  fluid  dynamics,  or  "hydrodynamic" 
motion.  The  first  logical  improvement  in  their  degree  of  sophistication  was  the 
inclusion  of  material  strength  in  shear.  This  was  performed  after  the  elastic- 
plastic  formulation  by  Wilkins  (Reference  2),  and  includes  the  von  Mises  stress- 
dependent  yield  criterion.  The  next  logical  step  in  upgrading  these  codes 
toward  more  realistic  descriptions  of  material  behavior  is  the  inclusion  of  a 
rate  dependent  description  of  dynamic  yielding. 

TWo  models  for  strain-rate  dependent  yielding  have  been  programmed  into  the 
PUFF  66  running  system,  and  either  or  both  may  be  employed  in  a  given  problem; 
however,  only  one  description  should  be  used  per  material.  The  first  model  is 
based  upon  the  formulation  proposed  by  Taylor  (Reference  3)  to  explain  dynamic 
yielding  through  a  phenomenological  description  of  dislocation  motion  in  metals. 
The  other  is  a  macroscopic  description  proposed  by  Alverson  and  Hannagud  (Refer¬ 
ence  4) .  The  adaptation  of  each  model  is  presented  in  the  following  sections 
with  comments  about  their  use  in  running  a  problem. 

The  models  were  incorporated  to  give  the  user  of  the  FUFF  66  codes  an  option 
that  would  allow  the  use  of  the  von  Mises  yield  criterion  or,  depending  upon 
one's  knowledge  of  the  material,  the  use  of  a  model  whose  degree  of  sophistica¬ 
tion  ranges  from  a  single  parameter  in  the  macroscopic  relaxation  model  to  a 
microscopic  model  based  upon  dislocation  dynamics.  The  changes  in  programming 
occur  in  the  stress-strain  calculation  of  the  HYDRO  subroutine;  this  is  shown 
by  the  flow  schematic  of  figure  1. 

The  following  sections  briefly  describe  the  theories  outlined  in  References 
3  and  4,  along  with  comments  about  the  effects  imposed  upon  the  wave  propagation 
calculations  by  these  descriptions.  These  effects  are  additionally  displayed  by 
the  results  of  selected  calculations  performed  with  each  model.  The  user  of 
these  material  descriptions  is  urged  to  read  References  3  and  4  before 
implementing  these  changes  to  the  PUFF  codes. 
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STRESS  STRAIN  CALCULATION 


ENTER 


CONTINUE 


Figure  1.  I  ow  Schematic  of  Stress  Strain  Calculation 
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SECTION  II 

DISLOCATION  DYNAMICS  MODEL 

The  dislocation  model  of  Gilman  and  Johnston  (Reference  5)  as  adapted  by 
Taylor  to  describe  dynamic  yielding  phenomena  is  a  more  sophisticated  treatment 
than  those  yield  descriptions  based  only  upon  stress  considerations.  This 
vicl  'ng  occurs  when  the  shear  stresses  on  linear  defects  in  a  material's  struc¬ 
ture  are  sufficient  to  cause  the  structure  to  deform  to  some  arrangement  that 
requires  a  lesser  stress  to  maintain  the  imposed  total  strain.  In  plastic 
deformation  of  crystalline  materials,  the  theory  of  dislocation  motion  was 
postulated,  and  has  since  been  verified,  as  the  mechanism  for  plastic  deformation 
of  materials  under  shear  stresses.  The  model  combines  elastic  response  with  a 
strain-rate  controlled  relaxation  to  hydrodynamic  or  plastic  behavior,  tyon 
initial  loading,  and  before  plastic  deformation  has  occurred,  the  material  has 
not  reconformed  its  structure  to  relieve  stress  loads,  and  the  shear  stresses 
are  at  their  maximum  value  for  the  imposed  total  strain.  The  average  disloca¬ 
tion  velocity  is  a  function  of  the  resolved  shear  stress,  and  the  rising  shear 
stress  produces  plastic  shear  strain  through  dislocation  motion.  The  defor.  j- 
tion  work  in  the  material  will  produce  new  dislocations,  and  the  resulting  higher 
dislocation  density  coupled  with  a  large  average  dislocation  velocity  produces 
the  rate  of  relaxation  of  the  total  stress.  When  the  relaxation  rate  and  the 
strain  rate  become  matched,  the  material  can  attain  an  equilibrium  stress-strain 
state.  This  state  will,  in  general,  lie  above  the  material's  hydrostat,  and  Is 
the  Hugoniot  state  one  would  measure  in  an  impact  experiment. 

For  this  model  Taylor  derives  the  relation  for  the  normal  stress  in  one- 
dimensional  strain  as 

o  -  (\  +  2u)  cT  -  |  u  y  (1) 

n  n  3 

where  o  is  the  normal  stress,  *  and  i.  are  the  Lamd  constraints,  is  the  total 
n  n 

normal  strain,  and  y  is  the  plastic  shear  strain.  By  using  the  expression  for 
the  time  dependence  of  the  plastic  shear  strain,  one  obtains 

-  bNv  (2) 
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where  b  is  the  Burgers  vector,  N  the  mobile  dislocation  density  or  the  total 
dislocation  length  per  cubic  centimeter  and  v  the  average  dislocation  velocity. 
The  relations  which  Taylor  uses  for  N  and  v  are 


N  +  ay 
o 


B 


(3) 


and 


tq  +  *y 


v  e 


(4) 


where  is  the  initial  dislocation  density,  a  and  6  are  constants  measured  for 
several  materials  by  Keh  and  Weissman  (Reference  6),  v  is  normally  taken  as  the 
acoustic  shear  velocity,  is  a  constant,  <f>  is  a  *  ork  hardening  term,  and  t  is 

the  resolved  shear  stress.  One  may  derive  expressions  for  the  plastic  shear 
strain  and  resolved  shear  stress  at  45°  to  the  propagation  direction  as 

Y  -  f^[(A+2u)e  -  oj  (5) 

and 

T  "  l[°  ~(X  +  f  u)e]  (6) 

Taylor's  equation  for  the  time  dependence  of  the  total  stress,  based  upon  stress 
being  positive  in  tension,  is 

0  "  -p  +  Sd 

as  in  the  Wilkins  formulation  for  the  elastic-plastic  flow  (Reference  2).  This 
ex^ressi  can  be  rewritten  to  allow  strain-rate  dependence  of  the  deviator 
only,  using  the  PUFF  66  convention  of  stress,  which  is  positive  in  compression. 
Since  the  hydrostatic  component  of  the  stress  may  be  written  as 

P  -  (a  +  |  p)  e 

the  stress  deviator  frrr.  Eq.  (1)  is  now  in  strain-rate  dependent  form: 


(7) 

a 


I 
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I 


» 


In  the  PUFF  formulation  t  becomes  -  —  and 

v 


By  combining  Eq.  (3)  and  (4)  with  the  expression  for  y  in  Eq.  (2),  one  has 


To+*Y 


By  using  the  PUFF  convention  for  calculating  stress  where 


•  9  » 

cf  *  P  —  S  , 

d 


this  equation  becomes 


To+1J>y 


By  incorporating  the  expression  for  y  and  t  In  Eq.  (5)  and  (6)  one  may  write 
the  equation  for  the  time-rate-of-change  of  the  deviatoric  stress  as 

v  !  ^  v +  f  ybv»jNo +  a[kt(x  +  2p)e"-J]  6jx  5 

S  «^exp  -|tq  +  |^[(A  +  2y)e-a]^ /~[o-(A  +  2/3  v)t]j  (8) 

1.  Mathematical  Dependence  of  Parameters 

The  physical  parameters  included  in  the  final  equation  for  this  model  are 
numerous,  and  one  might  expect  them  all  to  affect  the  calculations  in  a  manner 
that  would  be  difficult  to  allow  a  parametric  sensitivity  study  to  shed  light 
on  the  specific  effect  each  parameter  would  have  upon  the  calculations.  It  is 
a  good  idea,  therefore,  to  look  at  the  actual  expressions  in  some  detail  to  get 
a  better  understanding  of  the  effects  of  the  parameters. 

Starting  with  Eq.  (1),  one  sees  that  this  can  be  broken  down  to  a  thermo¬ 
dynamic  P-V-E  state  and  a  shear  strain-rate  dependent  deviatoric  stress  component 
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The  8/3  uy  term  simply  relates  the  extent  of  reconformation  of  material  structure 
to  the  rate  of  relaxation  of  the  total  normal  stress  in  a  one-dimensional  strain 
calculation.  Returning  to  the  initial  expression  for  shear  strain  rate, 

Y  *  bNv 

the  product  of  the  three  terms  is  simply  the  volume  of  slippage,  or  deformation, 
which  has  occurred  in  the  material  per  unit  volume,  per  unit  time.  One  must  now 
evaluate  the  shear  stress  and  shear  strain  dependence  of  N  and  v  to  determine 
the  shear  strain-rate. 

The  formulation  used  to  evaluate  the  average  dislocation  velocity  v,  is  given 

by  Eq.  (4).  This  is  a  continuously  varying  function  of  the  type  shown  in  figure 

2.  One  notes  that  the  average  dislocation  velocity  asymptotically  approaches 

the  value  v  ,  and  therefore  could  never  reach  this  maximum  in  a  normal  calcula- 
00 

tion.  The  work-hardening  term  in  Eq.  (4)  is  included  to  account  for  the  inter¬ 
actions  of  the  dislocations  with  each  other.  There  are  a  number  of  postulated 


Figure  2.  Average  Dislocation  Velocity  versus  Resolved  Shear  Stress 
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physical  mechanisms  that  may  occur;  however,  their  overall  effect  should  be 
thought  of  as  increasing  the  resolved  shear  stress  required  to  produce  a  given 
average  dislocation  velocity. 

The  change  in  the  number  of  dislocations  present  in  the  material  is  generally 
a  function  of  the  deformation  work  being  done  on  the  material,  and  is  expressed 
for  this  model  in  Eq.  (3)  as 

N  =  N  +  ay^ 

0 

The  coefficient  6  has  been  generally  determined  to  fall  between  the  values  of 
0.5  and  1.5  by  Keh  and  Weissman  (Reference  6).  As  the  plastic  shear  strain,  y, 
will  have  a  value  much  less  than  1,  one  can  see  that  the  rate  of  increase  in 
dislocation  density  is  inversely  proportional  to  3-  Thus  the  rate  of  relaxation 
is  inversely  proportional  to  3,  and  directly  proportional  to  a  if  all  other 
quantities  in  Eq.  (8)  are  fixed.  Since  this  is  not  the  case  in  a  dynamic 
calculation,  these  parameters,  along  with  the  others,  shou]d  only  be  evaluated 
in  a  qualitative  sense.  The  dependence  of  N  on  6  and  y  is  illustrated  in 
figure  3. 

One  should  notice  in  this  formulation  several  things  that  indicate  that  it 
is  restricted  to  use  in  compression,  pending  some  modifications.  First,  the 
plastic  shear  strain  should  be  written  as  a  differential  quantity  if  subsequent 
unloading  is  to  be  handled.  Second,  one  sees  that  the  stress  is  always  driven 
toward  tension  as  currently  formulated,  and  can  cause  the  unloading  to  be 
unbounded  without  constraints.  This  could  be  corrected  by  checking  to  determine 
if  the  material  is  unloading,  and  then  forcing  the  time  dependent  term  to  drive 
the  stress  deviator  toward  the  hydros  tat.  Another  method  would  be  to  use  the 
normal  von  Mlses  yield  criterion  in  unloading. 

2 .  Inclusion  of  Incremental  Shear  Strain 

By  allowing  the  calculation  to  handle  plastic  shear  strain  as  an  incrementally 
increasing  function,  one  eliminates  the  need  for  using  the  dislocation  model  only 
during  compressive  loading.  The  changes  necessary  to  perform  this  type  of 
calculation  are  relatively  straightforward  and  involve  only  a  few  additional 
statements  in  the  calculation  sequence.  The  terms  affected  are  the  plastic  shear 
strain  and  the  resolved  shear  stress  calculations. 
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The  resolved  shear  stress  as  currently  calculated  in  Eq.  (6)  must  be 
modified  to  become  positive  in  both  compressive  loading  and  unloading,  as  the 
expression  for  average  dislocation  velocity  otherwise  becomes  unbounded. 
Physically  this  modification  simply  says  that  the  amount  of  plastic  shear  strain 
is  a  scalar  quantity,  and  that  a  revolved  shear  stress,  tensile  or  compressive, 
results  in  dislocation  motion.  Equation  (6)  is  rewritten  as 


x  -  r  ABS 
4 


'(A  +f  u)e]or*  f  [Sd] 


The  plastic  shear  strain  is  now  a  quantity  that  must  increase  on  loading  and 
unloading,  and  is  calculated  by 


n+1  n  .  . 

Y  *  Y  +  Ay 


The  increase.  Ay,  must  therefore  be  calculated  at  each  time  increment,  n,  and 


can  be  written 


.  nH  *n+l  .  n+1 

Ay  «  y  •  At 


The  difficulty,  of  course,  is  in  determining  a  value  of  y  that  is  consistent 
with  the  values  of  Nn+3  and  vn+3.  This  could  be  an  iterative  process,  or  one 
could  assume  that  the  Ay  will  be  small  with  respect  to  y  and  proceed  in  the 
following  manner.  One  already  has  the  value  of  the  deviatoric  stress  without 
any  relaxation  from 

e  *  c  n  4  V  ..n+1 

Sd  Sd  +  3  U  v  At  (13) 


One  may  approximate  the  value  of  resolved  shear  stress  as 


n+1  3 

t  •  r  ABS 
4 


and  the  average  dislocation  velocity  as 


n+1  To  +  *Y 

v  -  v  exp - — t- 
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The  dislocation  density  is  obtained  by 

n”**  ■  nq  +  a(Yn) 

and  the  plastic  shear-strain-rate  by 

•n+1  ,  .Ji+l  n+1 

Y  ■  b  N  v 


(16) 


(17) 


From  this  one  obtains  the  "relaxed"  value  of  the  deviatoric  stress  and  the  new 
value  of  plastic  shear  strain  by  the  relations, 


and 


n+1 


+  f  u  Y 


n+1 


At 


n+1 


n+1 


n 


Y  +  y 


n+1  . kn+l 
At 


(18) 


3.  Comparison  with  Local*  Shear  Strain  Calculation 

One  may  wonder  if  this  calculation  and  the  one  -involving  the  local  shear 
strain  give  the  same  results.  To  make  this  determination,  two  test  problems 
with  identical  initial  conditions  and  material  properties  were  run.  Both  runs 
used  the  von  Mises  yield  criterion  in  relaxation,  making  the  two  problems  as 
identical  as  possible.  The  problem  used  a  l/8-inch  driver  plate  of  ARMCO  iron 
with  an  initial  velocity  of  1,55  x  104  cm/sec.  The  target  was  1/4-inch  of 
ARMCO  iron  backed  by  a  slab  of  quartz.  A  Lagrangian  edit  was  taken  just  inside 
the  quartz  at  the  iron/quartz  interface,  and  the  observed  waveform  was  used  for 
the  basis  of  comparison. 

The  material  properties  used  for  the  iron  were 

Bulk  modulus  •  9.05  x  1011  dynes/cm2 
Shear  modulus  ■  8.12  x  10 11  dynes /cm2 
Burgers  vector  ■  2.5  x  10”8  cm 
Shear  velocity  •  3.22  x  105  cm/sec 
Initial  dislocation  density  ■  106  lines/cm2 


♦Local  Shear  Strain  refers  to  that  calculated  at  any  instant  in  time,  Eq.  (5). 
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Coefficients : 

a  *  1.0  x  10 11  lines /cm2 
6  -  1.0 

tq  *  1.98  x  1010  dynes/cm2 
Initial  plastic  shear  strain  *  0. 

The  quartz  had  an  elastic  modulus  of  8.73  x  1011  dynes/cm2.  To  allow  the 
driver  to  separate  from  the  target,  an  interface  strength  of  5  dynes/cm2  was 
imposed. 

In  addition  the  incremental  plastic  shear  strain  model  was  employed  to 
control  the  deviatoric  stress  in  the  unloading  portion  of  the  stress  pulse.  The 
results  of  the  three  calculations  are  shown  in  figure  4. 

One  notes  some  differences  in  the  compressive  portion  of  the  wave  for  the 
local  and  incremental  shear  stress  models.  The  local  shear  stress  model  has  a 
higher  relaxation  rate,  causing  the  upper-lower  yield  phenomenon  to  be  less 
pronounced.  Because  the  Burgers  vector  is  constant,  and  v  is  calculated  in  the 
same  fashion  for  both  models,  the  only  parameter  that  can  be  controlling  the 
difference  is  the  value  of  the  plastic  shear  strain  in  the  relationship  for  the 
dislocation  density.  Thus  one  must  conclude  that  there  is  some  difference  in 
the  two  methods  used.  Because  the  difference  is  a  cumulative  effect  in  the 
Incremental  model,  the  magnitude  for  any  given  step  could  be  relatively  small 
and  difficult  to  resolve.  Also,  the  precision  of  the  models  used  does  not 
warrant  extensive  lnvestlgatior  of  these  differences. 

A  rather  significant  difference  is  noted  between  the  unloading  wave  profile 
for  the  von  Mlses  yield  criterion  and  the  rate  dependent  yield  function,  lhese 
differences  have  not  been  investigated  at  this  writing,  but  a  difference  would 
be  anticipated.  A  series  of  stress-versus-dlstance  plots  for  different  times 
are  shown  in  the  Appendix  for  the  total  plastic  shear  strain  model  with  the  von 
Mlses  yield  criterion  in  unloading,  and  the  differential  plastic  shear  strain 
model  for  both  loading  and  unloading.  This  gives  the  reader  a  better  example 
of  the  total  problem  that  produced  the  waveforms  in  figure  4. 
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Figure  4.  Comparison  Between  Waveforms  Generated  bv  Various  Dislocation  Models 
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SECTION  III 

STRESS-RELAXATION  MODEL 

Two  models  were  developed  in  reference  3.  One  was  a  macroscopic  model 
based  upon  a  single  characteristic  relaxation  time;  the  other,  a  microscopic 
model,  used  dislocation  dynamics  coupled  with  fairly  complicated  formulas  for 
the  rate  of  mobile  and  total  dislocation  production  and  disappearance.  While 
the  latter  model  is  somewhat  different  from  Taylor  s,  they  both  obtain  y  from 
the  mobile  dislocation  density  and  the  average  dislocation  velocity.  Thus,  the 
other  model  chosen  for  programing  into  PUFF  was  the  macroscopic  relaxation  model 
that  employs  a  single  parameter,  the  characteristic  relaxation  time,  to  produce 
a  time-dependent  stress  deviator  calculation. 

This  model  uses  the  difference  in  the  actual  devlatoric  stress  computed 
and  the  equilibrium  value  calculated  from  the  von  Mises  yield  criterion  along 
with  the  characteristic  stress  relaxation  time.  This  formulation  is  simply 
stated  as 


Act 

T 


(19) 


where  Ao  is  the  difference  between  the  devlatoric  stress  and  the  equilibrium 
stress  deviator  baaed  on  a  stress-dependent  yield  criterion,  and  T  is  the 
characteristic  relaxation  time.  One  must  take  into  account  whether  or  not  the 
devlatoric  stress  is  above  the  equilibrium  yield  value,  and  whether  or  not  the 
material  is  loading  or  unloading.  This  is  controlled  by  first  checking  the 
stress  deviator  in  the  normal  manner  to  determine  whether  the  yield  condition  is 
exceeded,  and  then,  if  the  time-dependent  relaxation  is  to  be  used,  determining 
the  new  value  of  the  deviator  based  upon  Eq.  (19).  This  equation  is  written  for 
a  given  lagrangian  zone  in  PUFF  as 

c  n+1/2  c  n  4  v  n+1 
S.  *  S  ,  +  r  ii  -  it 
d  d  3  v 


n.l  "*1/2  /  n+1/2  _  2  /t 

d  5d  \5d  3  7  T 


(20) 
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The  variable  0  is  +1  on  compression  and  -1  on  unloading.  This  is  determined 
within  the  stress-strain  calculation,  and  is  positive  when  the  total  stress  is 
above  the  hydrostat.  An  interesting  facet  of  this  model  is  that  it  can  be 
broken  into  a  hydrostatic,  elastic  deviatoric,  and  plastic  shear  contribution 
to  the  total  stress.  One  remembers  from  Eq.  (1)  that  the  time  dependence  of 
the  total  longitudinal  stress  was  written  as 


•  /  2  \  •  4  •  8 

°  "\X  +  3  u)e  +  3  u  e  "  3  u  Y 


whereas  in  this  model  it  becomes 


•  /  2  \  •  4 

a-^X+-jujc+-ue 


(V  \  y) 


One  may  relate  the  last  terms  in  Eq.  (1)  and  (21)  to  find: 


a  •  ( 

-  II  V  ■  - 


S  .  -  0 
a 


or  that 


;  .  i_p 

y  8w  L 


SD  -  9  2/3 
T 


Y  •  o  for 


d|  -  3 


2 

1 7  y 


The  plastic  shear  strain  is  given  in  Eq.  (5)  as 


i  (' +  2u)  ■  - '] 


so  that  the  plastic  shear  strain  rate  is 


’  *  at  ('  *  2l)‘  -  [ 


One  can  see  that  there  is  dimensional  equivalence  between  equations  23  and  24; 
however,  the  time  dependence  of  the  quantities  is  based  upon  different  physical 
properties.  One  can  also  envision  a  value  of  T  which  causes  the  relaxation  rate 
to  become  too  large  for  the  stability  criterion  in  PUFF.  This  will  occur  when 
T  is  smaller  than  the  time  step  in  the  problem.  Thus  one  must  use  caution  in 
choosing  T,  to  insure  that  it  is  always  greater  than  the  antic'pated  tine  step. 


14 


AFWL-TR-69-11 


i 


1 .  Test  Problems 

The  same  impact  conditions  were  used  here  with  the  comparison  runs  as  were 
performed  on  the  dislocation  model.  The  relaxation  time  was  varied  between 
10“4  and  10“8  seconds,  and  the  variation  in  the  waveshapes  measured  at  Lagran- 
gian  edit  are  shown  in  figure  5.  Since  one  must  keep  the  relaxation  time  larger 
than  the  time  step  for  stability,  no  lower  values  were  used.  One  observes  that 
these  effects  are  somewhat  different  than  the  effects  seen  from  the  dislocation 
model . 

2 .  Programming  Changes  in  PUFF 

Relatively  few  changes  are  required  in  the  PUFF  code  to  input  the  two  models 
discussed.  The  material  variables  must  be  read  in  with  the  GENRAT  routine,  and 
placed  in  the  common  block.  The  only  changes  required  in  the  program  are  in 
the  HYDRO  subroutine,  and  occur  in  the  stress-strain  calculation.  The  parameters 
used  are 

BURG(M)  ■  Burgers  vector,  (lines/cm2) 

CSHR(M)  ■  Acoustic  shear  velocity,  (cm/sec) 

TAUO(M)  ■  tq  in  equation  4,  (dynes /cm2) 

ALFR(M)  ■  a  in  equation  3,  (lines/cm2) 

BETA(M)  ■  8  in  equation  3,  (dimensionless) 

ANO(M)  *  Initial  dislocation  density,  (lines/cm2) 

FEE(M)  •  $  in  equation  4,  (dynes/cm2) 

TRELAX(M)  «■  Characteristic  relaxation  time  for  macroscopic  model,  (sec) 

The  changes  to  the  stress-strain  calculation  are  shown  in  figure  6;  the  general 
flow  of  the  calculation  is  shown  in  figure  1.  The  new  variable  names  were 
chosen  for  ease  of  understanding  with  respect  to  the  equations  previously 
discussed . 
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SECTION  IV 


DISCUSSION 

The  models  incorporated  into  the  code  may  serve  as  a  skeletal  framework 
that  may  be  varied  to  approximate  different  phenomena  relating  to  the  strain 
rate  dependent  behavior  of  materials  under  impulsive  loads.  The  macroscopic 
model  may  be  employed  either  as  used  for  the  runs  shown  in  this  report  with  the 
yield  strength  below  the  hydrostat  set  to  zero,  or  for  any  other  variation 
desired.  One  may  also  choose  to  make  the  yield  strength  time  dependent  as 
discussed  in  reference  3.  These  changes  are  easily  made  in  the  stress-strain 
calculation  presented  in  figure  6. 

The  dislocation  dynamics  model  is  based  upon  Eq.  (2),  and  the  expressions 
used  for  determining  the  dislocation  density  (mobile  or  total)  and  the  disloca¬ 
tion  velocity  (mobile  or  average)  may  be  modified  by  changing  the  expressions 
for  the  variables  ANI  and  VDIS,  respectively.  The  initial  value  of  the  plastic 
shear  strain  and  the  dislocation  density  should  certainly  be  properly  matched 
if  they  are  to  be  used  meaningfully.  One  can  get  some  idea  of  the  effect  upon 
the  calculation  if  the  initial  plastic  shear  strain  is  non-zero  by  examining 
figure  7.  Here  a  comparison  is  made  between  the  Lagrangian  edit  in  the  quartz 
for  runs  using  initial  plastic  shear  strain  values  of  zero  and  1  percent.  Bv 
allowing  the  material  to  have  been  previously  "cold  worked"  to  the  1  percent 
value,  a  common  practice  with  mild  steels,  one  finds  that  the  upper-lower 
yield  phenomenon  no  longer  appears  in  the  propagated  stress  wave,  The  effects 
of  the  various  other  parameters  upon  a  waveshape  have  been  lightly  touched  on 
previously;  however,  the  user  should  familiarize  himself  with  their  significance 
by  performing  a  comprehensive  sensitivity  study  on  the  parameters  of  whichever 
model  he  chooses  to  employ. 


-  INITIAL  PLASTIC  SHEAR  STRAIN  =  0  PERCENT 
-INITIAL  PLASTIC  SHEAR  STRAIN=.OI  PERCENT 
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Figure  7.  Illustration  of  the  Effect  of  the  Initial  Value 
of  Plastic  Shear  Strain  on  a  Propagated  Waveform 
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SECTION  V 
CONCLUSIONS 


Based  upon  the  computed  stress  waves  presented  for  the  two  strain-rate 
dependent  material  response  models  employed,  the  following  conclusions  are 
drawn : 

1.  More  realistic  descriptions  of  the  response  of  crystalline  materials  are 
possible  by  employing  deviatoric  stress  calculations,  to  account  for  microscopic 
behavior  in  the  PUFF  code. 

2.  The  flexibility  of  the  models  presented  provides  a  more  physically  realistic 
description  of  material  response  of  impulsive  loads  and  necessitates  more  precise 
definition  of  the  material's  microscopic  behavior  for  these  descriptions  to  have 
any  real  meaning. 

3.  The  microscopic  model  employs  enough  parameters  to  make  concise  definition 
of  the  effect  of  any  single  parameter  upon  a  propagated  stress  wave  virtually 
impossible. 
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APPENDIX 

Two  of  the  computer  runs  described  in  the  report  are  more  fully  presented 
here.  A  series  of  stress-versus-distance  waveforms  are  displayed  for  both  the 
incremental  plastic  shear  strain  model  used  for  loading  and  unloading,  and  tor 
the  local  plastic  shear  strain  model.  The  cycle  numbers  and  problem  times 
associated  with  each  plot  appear  at  the  top  of  each  margin.  Figures  8  through 

12  are  the  results  for  the  incremental  plastic  shear  strain  problem,  and  figures 

13  through  17  are  from  the  local  plastic  shear  strain  problems. 
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(B)  CYCLE  1175  TIME  =  1.686  MICSEC 


Figure  10.  Stress-Distance  Histories  for  Incremental  Shear  Strain  Proble 
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